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Abstract 

The noiseless limit of a thermal ratchet device using ferrofluids is studied in detail. 
Contrary to previous claims it is proved that no directed transport can occur in this 
model in the absence of fluctuations. 
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1 Introduction 



Rectification of non-equilibrium fiuctuations can be accomplished with the 
help of so-called ratchets [1]. In these devices a periodic potential, i.e. force 
field with zero spatial average, and undirected random noise conspire to pro- 
duce directed transport. Besides their fundamental importance for statistical 
mechanics [2] ratchets have gained renewed interest under the name of "Brown- 
ian motors" due to their possible relevance for transport in biological cells and 
potential applications in the field of nano-technology [3]. For a comprehensive 
review of the field see [4]. 

Recently a thermal ratchet system using ferrofiuids was introduced [5,6]. Fer- 
rofiuids are colloidal suspensions of ferromagnetic grains in a suitable carrier 
liquid [7]. The spatial orientation of the ferromagnetic particles is influenced 
by the local vorticity of the flow field of the carrier liquid as well as by thermal 
fiuctuations due to random collisions with the molecules of the liquid [8]. More- 
over this orientation can be coupled via the magnetic moment of the grains 
to external magnetic fields. Choosing a suitable time dependence of this field 
to drive the system away from equilibrium it is possible to rectify the orienta- 
tional fiuctuations of the ferromagnetic particles. More precisely an external 
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magnetic field without net rotating component can be used to set up noise- 
induced rotations of the ferromagnetic grains. Besides other advantages this 
system has the attractive feature that by the viscous coupling of the particles 
to the surrounding liquid the angular momentum of the many nanoscopic mo- 
tors is transferred to the carrier liquid and shows up as a macroscopic torque 
density which can be easily detected in experiments [5]. 

The necessity of thermal fluctuations for the operation of the ratchet of the 
described type in ferrofluids was disputed in [9]. In fact there are several 
examples, e.g. so-called rocking ratchets (see [4]), in which for appropriate 
choices of parameters directed transport may even occur without fluctuations. 
However, in order to spin up ferrofluid particles by a magnetic field without net 
rotating component as introduced in [5] the presence of thermal fluctuations is 
indeed indispensable. This is shown in the present paper where we prove that 
in the deterministic dynamics no full rotations of the particles may occur and 
that no torque may be transferred from the magnetic field to the particles. 

To this end we first recall in section II the basic equations for a single ferro- 
magnetic particle in an oscillating external field as derived in [5]. The ratchet 
effect was shown to occur in strongly diluted ferrofluids also such that dipole- 
dipole interactions between the ferromagnetic grains may safely be neglected 
and a single-particle picture is appropriate for the theoretical analysis. Setting 
the noise intensity equal to zero we investigate in sections III and IV the de- 
tails of the deterministic dynamics of the particle and show that only solutions 
without full rotation of the particle are possible. Finally, section V contains 
the main conclusion. 



2 Basic equations 

We consider a spherical particle of volume V and magnetic moment m, subject 
to a time dependent magnetic field of the form 



where Hy{t) is a periodic function with period 2tt/uj. As instructive example 
[5] we will consider the special case 



where a, /? and the phase shift 6 are control parameters. The particle is im- 
mersed in a fluid of viscosity 77. 

To describe the orientation of the particle we use the unit vector e = m/m 
where m denotes the the magnetic moment of the particle and m its modulus. 



(1) 



Hy{t) = a cos{ujt) + p sin(2u;t + S) 



(2) 
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Changes of e are described by the equation 



-=nxe. (3) 
where is the angular velocity of the particle. 

Furthermore we consider an overdamped stochastic dynamics in which the 
magnetic torque 

^mag = me X H (4) 

and the stochastic torque [11], which results from the interaction between the 
particle and the surrounding liquid, 

Nstoch = V2D^it) (5) 

is counterbalanced by the viscous torque [10]: 

^vis = -6vvn (6) 

In equation (5), ^{t) is a vector of Gaussian white noise with zero mean and 
unit variance. The noise intensity D is related to the temperature T of the liq- 
uid by the Einstein relation: D = 6r]VkBT, where /c^ stands for the Boltzmann 
constant. From (4), (5) and (6) we find: 

6r]Vn = mexB. + V2D^{t). (7) 

This relation together with equation (3) yields a closed equation for the time 
evolution of e of the form 

de m , \/2-D ^ , , 

— = (exH)xe + - £ x e. (8) 

Introducing dimensionless units we measure time in units of the inverse driv- 
ing frequency, t t/u^ and use QrjVuj/m as unit for the magnetic field 
strength H — » {QrjVu /m)'iii. The noise intensity D is scaled according to 
D {QrjV)^ D. Eq. (8) then reduces to 

— = (e X H) X e + V2D £ x e. (9) 

(JjV 

It is convenient to parametrize the orientation of the particle by the two angles 
6 and according to 

e = (sin ^ cos 0, sin 6* sine/), cos ^) (10) 
The Langevin equations for the time evolution of these angles are then of the 
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form [11,12] 



de 

dt 

d(j) 

It 



~U + Dcote + V2D^e 



d 



sin" e I 



u 



sin 9 



(11) 
(12) 



where we have introduced the potential, 



U{d,(t)) = -e ■ H = - sin 6* (if^ cos </) + Hy{t) sin( 



(13) 



The observable of principal interest for the thermal ratchet effect in ferrofluids 
is the average torque N arising at the magnetic particle in the long time limit. 
Here the average is over time and hence includes both the ensemble average 
over different realizations of the noise as well as the average over the time 
dependence of the external magnetic field. The focus of the present paper is 
on the T ^ limit implying D — ^> 0. The other system parameters like the 
fluid viscosity are assumed to stay constant. We then find for the averaged 
torque from the dimensionless forms of (4) and (6): 



N = exH = -ri, (14) 

since the stochastic torque is zero in the absence of fluctuations. 

For later use it is instructive also to study the case in which the particle 
orientation is confined to the plane defined by 6 = 7i/2. In this case we have 
N = (0,0, iVj and Q = (O,O,90/9t). Therefore from eq. (14) we find 

N: = lim f'dt' ^{t') (15) 

and hence 

(t2— tl)— >oo t2 — ti 

In the absence of particle rotation, i.e. if < < 27r, therefore no 
torque may arise. 

In the following two sections we study the deterministic dynamics given by 
(9) with = in detail. We first consider the case = and then deal with 
the more general situation where 0. 
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3 The case H,, = 



In this case the external field has the form 

H = (0,/J,(t),0) (17) 

with a general time-dependent Hy(t). For D = eq. (9) yields the following 
equations for the components of e: 

de 

= -e^Hy{t)ey (18) 

dp 

= -elHy{t) + Hy{t) (19) 

de 

= -e,Hy{t)ey . (20) 

(21) 

This system of differential equations is to be completed by appropriate initial 
conditions at some initial time to. Without loss of generality we can choose 
the coordinate system in such a way that 62(^0) = 0, i.e. we take the x-2;-plane 
as the plane defined by e(to) and the direction of the magnetic field. From 
equation (20) then follows that is identically zero, Czit) = 0. 

>From (18) and (19) we find 

ey{t) = tanh (^J^ dt'Hy{t') + arctanh ej/(to)) (22) 

and ^ 

ex(t) = e,(to) exp ^ dt' Hy{t')ey{t')^ (23) 

The integral / in the exponential function in (23) can be determined using 
(22) 

I = dt'Hyit') tanh dt"Hy{t") + arctanh Cyito)^ . (24) 

Substituting u{t) = //^ dt'Hy{t') + arctanh ey(to) this gives 

r^*^ , w N , / cosh(n(t)) \ 
/= / du tanhiu) = In \ — , / ; 25 

Ju(to) ^ ' \cosh{u{to))) ^ ' 

and we finally get the solution 



cosh ( arctanhe„(to 

ex{t) = e^{to) 7 ^ -—^ ^ (26) 

cosh f arctanh (to) + Jto dt'Hy{t') 



ey(t) = tanh (it'i7y(t') + arctanh ey(to)^ • (27) 



5 



-0.5 I 1 1 1 1 " 

5 10 ^ 15 20 25 

Fig. 1. Trajectories (p{t) for the case ffx- = and for the time dependence of 
Hy{t) as defined in (2) and to = 0. The parameter values are chosen as a = 1, 
(5 = 1, and 5 = 0. The different curves correspond to different initial conditions, 
4>o = 7r,57r/6,7r/6 and (from top to bottom). There are no rotary solutions in the 
absence of fluctuations and hence no average torque may arise. Note that the trajec- 
tories are different for different initial conditions, there is hence no unique long-time 
behaviour. 



>From (26) it follov^^s that there are no full rotations of the particle, because 
the x-component of e cannot change sign. Expressed in terms of and 0, 
equations (26), (27) for this choice of coordinates become 




/sinh (jl dt'Hy{t' 
(f)(t) = arctan 



where 0o = 0(^o) The domain of the arctan function has to be chosen such 
that arctan (tan(0o)) = 0o- Some trajectories of 0(t), for the time dependence 
Hy(t) as given in (2) and a certain choice of parameters are shown in figure 1. 
Note that there is no unique long time behaviour for = 0. The trajectories 
of (j){t) are different for different initial conditions. 

The main conclusion of this section is the absence of full rotations of the 
ferrofluid particle for = 0. Hence is bounded to an interval 0mm < 
< <Pmax- Together with eq. (16) it follows that for = and any time 
dependence Hy{t) there is no average torque. 
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4 The case H,, > 



In the previous section we have discussed the behavior of a ferromagnetic 
particle in a time-dependent field in ^/-direction and we have shown that in 
this case no full rotations of the particle may occur. Intuitively it is clear 
that an additional constant field in x-direction should not change this result. 
Nevertheless the argument given in the previous section according to which 
the a;-component of the orientation vector e cannot change sign does no longer 
hold true. If, e.g., e^ito) is negative and is positive while Hy{t) = 0, the 
orientation of e will tend toward the direction of H and therefore has to 
change sign. It is therefore necessary to take a closer look on the possible 
implications of a constant magnetic field in x-direction for the motion of the 
magnetic particle. We will show in this section that even in the presence of 
such a field 6 will converge to 7r/2 in the long time limit while will be confined 
to an interval (pmin < < (pmax- 

We start with eqs. (11) and (12) in the deterministic limit, D = 0. Using 
eq. (13) they can be written in the form 

d6 

— = cos 6 (^Hx cos 4> + Hy{t) sin (30) 

^ = -^H, (Git) cos - sin 0) , (31) 
at smd ^ ' 

where we have introduced the function Git) = Hy(t)/H.j.. 
Let us first look at the sign of dcp/dt which is given by 

sgn iy^^ = sgn(cos0) sgn(G'(t) — tan0). (32) 

Since G{t) is a periodic function it has a maximum Gmax and a minimum 
Grain- Deuotc by (pmax and 0mm the solutions of the equation tan0 = Gmax and 
tan0 = Grain respectively in the interval (— 7r/2, 7r/2). Assuming to belong 
to the interval (— 7r/2, 37r/2) it is easy then to show that d(j)/dt is always (i.e. 
for all t) positive if —7i/2 < < (prnin or (prnax + vr < < 37r/2 (region II), and 
that it is always negative for (prnax < < 0mm + 7r (region I) (cf. fig. 2). For the 
remaining values of the sign of d(j)/dt depends on the actual value of G{t). 
Note that at = ±7r/2 both {G{t) — tan0) and cos0 change sign such that 
dcp/dt does not. Hence these "critical" points corresponding to e^; = belong 
to regions in which the sign of d(f)/dt is independent of time. 

In order to discuss now the time evolution of 9 and we have to consider 
different initial conditions for 0. Quite generally we may assume H^c > 
without loss of generality. Let us first consider the case —7i/2 < 0o < 7r/2, 
i.e. starts in regions I, II, or III. Then 0(t) has to reach region III sooner 
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Fig. 2. Left: Regions with different values of dcp/dt as shown by the arrows. In region 
I d(j)/dt < for all t, in region II d(f)/dt > for all t. In regions III and IV the sign 
of d(j)/dt depends on the actual value of G{t). Right: Definition of (l)max and 4>min 
and graphical determination of the sign of d(j)/dt. 

or later and will not be able to leave it again (cf. fig. 2). Strictly speaking this 
is correct only if (f){t) evolves continuously. But we have to take into account 
also that there is the possibility of discontinuous changes of when 6 reaches 
the limiting values 9 = or 9 = n. However, this cannot happen either. To 
prove this statement we introduce the quantity 



h = -. 

Using (9) with = its time derivative is given by 



dh _ 
lit ~ 

or expressed in terms of 9 and <; 

dh 
Itt 



Hx 1 
-h. 



cos d) sin 9 



h. 



(33) 



(34) 



(35) 



Hence 



h{t) = h(to) exp 



- / dt' 

Jto cos0(t')sin^(f) 



(36) 



Accordingly, as long as stays in the interval (— 7r/2..7r/2) the integral in the 
exponent of (36) grows and consequently \h(t)\ monotonically decreases with 
time. Therefore 9 cannot reach or vr. 



Hence (j){t) cannot leave region HI neither by continuous nor by discontinuous 
changes. On the one hand this ensures that (pmin < (p(t) < 'Pmax for all t on 
the other hand it implies via (36) that h{t) for large t and hence that 9 
converges asymptotically to 7r/2. 
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If 7r/2 < 00 < 37r/2 the evolution of starts in regions I, II, or IV. If at some 
later time 0(t) is found in the interval (— 7r/2, 7r/2) we are back to the previous 
case. If not and /i(to) 7^ eq. (36) implies that \h{t)\ increases monotonically 
with time and accordingly 6 tends to either or vr. By symmetry both cases are 
equivalent so each other so let us focus on 6* ^ 0. Then h 1/6 and eq. (35) 
acquires the asymptotic form dh/dt = C h"^ with some positive constant C. 
Therefore there will be a finite time singularity in the solution h{t) and we 
get 0{ti) = for some finite ti. Since the magnetic field has a positive x- 
component it is clear that for t > ti we will have —tt/2 < (f){t) < 7r/2 and 
hence we are again back to the first case. 

Summing up, except for a set of measure zero, namely 0mm +7r < 0o < (pmax+T^ 
and 9 = 7r/2, all initial conditions give rise to a long time dynamics with values 
of between 0mm and (pmax- In any case also for 7^ no full rotations of 
the particle are possible since these would imply that 0(t) lies for some t in 
region III which it were unable to leave again. As in the case = we then 
find from (16) that no average torque is transferred from the magnetic field 
to the particle or its surrounding liquid. 



5 Conclusion 

By a detailed analysis of the deterministic dynamics of a magnetic dipole in 
an external magnetic field with constant x-component and time periodic y- 
component we have proved that no full rotations of the particle may occur. 
This shows that for the thermal ratchet effect in ferrofluids reported in [5] 
thermal fluctuations are indispensable. 
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